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Abstract 

We study a mean value of the classical additive divisor problem, that 
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d{n + l)d(n + / + /)— main term , 

with quantities iV > 1, 1 < F < iV^"= and 1 < L < iV. The main term 
we are interested in here is the one by Motohashi [27] , but we also give 
an upper bound for the case where the main term is that of Atkinson 
pp. Furthermore, we point out that the proof yields an analogous upper 
bound for a shifted convolution sum over Fourier coefficients of a fixed 
holomorphic cusp form in mean. 
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1 Introduction 

A widely studied occurrence of the classical divisor function d{n) — X]rf|n 1 
the additive divisor problem, in which one investigates the asymptotic behavior 
of the sum 

D(x-, m) — d{n)d{n + m) (1) 

as X tends to infinity and to is a given positive integer. Ingham was the first to 
find an asymptotic formula for the above snm in 1927, and in 1931 Estermann 
brought the Kloosterman sums into discussion, improving the earlier result. The 
importance of uniformity in the shift m was noted ten years later by Atkinson 
[T], who further used his findings in studying the fourth power mean of the 
Riemann zeta-function on the critical line. In 1982 Deshouillers and Iwaniec 
applied Kuznetsov's trace formulas into the problem, obtaining a notable im- 
provement to the upper bound of the error term. Finally, especially crucial for 
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this paper, Motohashi [57] worked out an explicit spectral decomposition for 
a weighted convolution sum, separating a main term and obtaining an upper 
bound for the error term uniformly in the shift m. We mention shortly papers 
by Heath-Brown, Iwaniec |TI] and Ivic- Motohashi [T^, who obtained mean 
value results for the error term, and and [52] for other related results. 

For a more thorough discussion of the history and the references of the above 
mentioned classical papers we refer to Motohashi's comprehensive paper |27j . 

In recent years the analogy between the divisor function and the Fourier 
coefficients of holomorphic and non-holomorphic cusp forms has appeared re- 
peatedly in the literature. The behaviour of the analogous sums to ([T]), 

A{x; m) = a(n)a{n + m) (2) 

wCx 

involving the Fourier coefficients of a holomorphic cusp form and 

= ^ t{n)t{n + to) (3) 

involving the Hecke eigenvalues corresponding to Fourier coefficients of a Maass 
form, has been studied intensively. In 1965 Selberg studied the meromorphity of 
a certain Dirichlet series related to the Fourier coefficients of holomorphic cusp 
forms, while Good in 1981 estimated the second moment of a modular L-function 
on the critical line using the spectral decomposition of shifted convolution sums 
coming from the Fourier coefficients of a holomorphic cusp form. Jutila |14j . 
[15] derived explicit asymptotic formulae for all three cases ([1]), ([2]) and ([3]) in a 
unified way, using the respective generating Dirichlet series. For a discussion of 
the substantial history including e.g. the works by Blomer, Duke, Friedlander, 
Harcos, Iwaniec, Jutila, Michel, Motohashi and Sarnak as well as a general 
spectral decomposition for the weighted shifted convolution sums we refer to 
the paper [2] by Blomer and Harcos. Along with these references we mention 
also 0, [7], [T6j and [2QJ. Together the three sums above go under the name 
of the shifted convolution problem. 

As a prologue for this paper, in Lemma 3 of his paper [16 Jutila studies the 
sum ^ in mean, and his argument can easily be extended to prove 
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0</<Fl<n<Af 1<1<L 
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a{n + l)a{n + f + I) 



<^ {N + F)'' N'' L (4) 



for all N,F > 1 and 1 < L < N . Instantaneously the proof gives also an upper 
bound for a sum over Hecke eigenvalues of Maass forms in mean, as stated in 
Lemma 6 of 03] : For iV, F > 1 and 1 < L < iV 



E E E tin + l)t{n + f + l) 

0<f<F l<n<N 
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1+s j 



<t: {N + Fy+^NL. (5) 



The idea in the estimation of the triple sum lies in its sensitiveness of the size 
of the innermost sum over I. A similar result, sensitive for the size of the shift, 
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is needed in the doctoral thesis of the author ([3T], Lemma 3.4). An analogous 
result over more general settings also appears on p. 81 in the paper [4] by 
Blomer, Harcos and Michel, leading to the same bounds (|1]), ([S]). In all the 
three papers [T^, [T7], [1] the motivation behind these studies has been gaining 
information about the upper bounds of the original shifted convolution sums 
and ([3]), while the author needed her estimate in her thesis in estimating a 
certain spectral sum over inner products involving holomorphic cusp forms and 
Maass forms. 

However, to our knowledge there does not exist a similar bound for the case 
of the divisor function, and the extension of the proofs of the cuspidal cases seem 
problematic due to the main term of the additive divisor problem. This paper 
is our attempt to bring new light for this case. Hence we state the following 
Theorem: 

Theorem 1.1 (Main theorem). Let N > 1, I < L < N and 1 < F ^ N'^'" . 
Then 



d{n + l)d{n + l + f) - — / 
Ir^L J(L+n)/f 



{2L+n)/f 



m{x] f)dx 



N^+^LF. 



(6) 



Here m(x; f) is as in (1.12) in J^Tj/ . 

The notation m ^ M stands for M < m < 2M . For the classical main term 
M{N, /) by Atkinson ([IJ, p. 185) we easily obtain the following upper bound: 

Corollary 1.2. Let N > 1, I < L < N and I < F ^ N^-^ . Then 
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fr^Fnr^N 



d{n + l)d{n + l + f)- M{2L + nj) + M{L + n, f) 



< Ar2+e + N^+^LF + N-'^+^L'^F^. 



Remark 1.3. With small values of the quantities F and that is, if L^F = 
o{N), the trivial bound N^'^'^L^F yields a better result than Theorem 11.11 and 
Corollarv ll.2l On the other hand, in their paper ^Oj Ivic and Motohashi gain 
the upper bound 
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Yd{n)d{n + f)-- / 



N/f 



m{x] f)dx 



< ]Sf4:/3+epl/3 



uniformly for F < N^^^ ^ . Here, according to Meurman ([25j, pp. 225-226) it 
seems possible to extend the range to _F < iV by replacing the upper bound by 
jy4/3+e^i/3 j^p^ rpj^j^ ^^^^j^ yields, up to factor 7V^ the same upper bound 
as in Theorem 11.11 in case F ^ N^^^ and L — N. Other than this, to our 
knowledge, our result ([6]) and the trivial bound together yield the genuinely 
best upper bound with the restriction F ^ A^^^^. 
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The main novelty in this paper is the methodological uniformity allowing an 
analogous treatment for the cuspidal cases: Following the analogous proof for 
the holomorphic cusp forms we immediately gain 



Theorem 1.4. Let N > I, I < L < N and 1 < F <. N^-" . The 
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To compare our result to the earlier mentioned, we notice that up to a factor 
and with the restriction F <ti N^~'^, ^ yields our result only if L x 1, and 
otherwise our result is genuinely better. 

In the case of non-holomorphic cusp forms we face the problem of the lack of 
a proper analogy for the spectral decomposition of the shifted convolution sum 
in question. There exists an analogous result for our crucial Lemma 12.31 namely 
Lemma 5 in [52], originally due to Jutila [13], [IH]. However, although the 
analogous part of this decomposition results the anticipated analogous bound, 
there appears also an arithmetic correction term, which leads to additional 
problems. If we were dealing with a genuinely oscillating weight function in our 
proof, then the contribution of the correction term would easily be proven to 
stay under the expected bound. However, as this is not the case we shall lighten 
this paper and just assign the conjecture below, leaving the proof for future 
study. 



Conjecture 1.5. Let N >1, 1 < L < N and 1 < F <s: 



Then 



EE 



J2t{n + l)t{n + l + f) 



N^+^LF. 



Remark 1.6. Note also Theorem 1 in [2], which gives a more general spectral 
decomposition of a shifted convolution sum including Hecke eigenvalues of the 
underlying cusp forms of two arbitrary cuspidal automorphic representations of 
PGL2(Z)\PGX2(K)- 



2 The needed notation and auxiliary lemmas 
2.1 The cusp forms 

Besides their analogous behavior to the divisor function, the holomorphic and 
non-holomorphic cusp forms play a role in the spectral decomposition of the 
shifted convolution sum over the divisor function (Lemma 12. 3p . Hence we start 
by introducing briefly some results concerning them. For the proofs and for a 
general reference the reader is referred to Motohashi's monograph |29j . 

We may confine ourselves to the cusp forms for the full modular group F = 
SL2 (Z) operating through Mobius transformations on the upper half plane H. 
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A holomorphic cusp form F{z) : H — > C of weight k E Z, with respect to F can 
be represented by its Fourier series 

oo 

F{z) = a{n)e{nz), e{a) = exp(27ria). 

n=l 

We may assume that k is even and k > 12, otherwise F{z) is trivial. We let 

{^j,k 1 1 < J < 

be an orthonormal basis of the unitary space of holomorphic cusp forms of 
weight fc, and write 

oo 
n=l 

We may suppose that the basis vectors are eigenfunctions of the Hecke operators 
Tk{n) for all positive integers n. Thus, in particular, Tk{n)ipj^k = tj,k{n)'4'j,k 
for certain real numbers tj_k{n), which we call Hecke eigenvalues. Comparing 
Fourier coefficients on both sides, one may verify that pj^k{n) = Pj^k{^)'tj,k{Ti) 
for all n > 1, 1 < j < d{k). We put a^-fe = 16F(2fc)(47r)-2fc-i x|pj-fc(l)|2, 
whence 

§(k) 

aj,k < k. (7) 

Furthermore, we let 

oo 

H],k{s) = y^Jj,kin)n~'' 

n=l 

stand for the Hecke L-function attached to ipj,k- The series converges absolutely 
for Re s > 1 because of the bound 

\tj,k{ri)\ < din) <^ (8) 

by Deligne [S]. Furthermore Hj^k{s) can be continued to an entire function, and 
it satisfies a functional equation, which implies that for bounded s 

H,,k{s) « (9) 

uniformly in j. Here c is a suitable constant, which depends only on Re s. 

A non-holomorphic cusp form u{z) = u{x + iy) : H — >■ C is a non-constant 
real-analytic F-invariant function in the upper half-plane, square- integrable with 
respect to the hyperbolic measure dfi{z) = over a fundamental domain of 
F. Furthermore u(z) is an eigenfunction of the non-euclidean Laplacian A = 
—y'^{-^ + ^ii"^ the corresponding eigenvalue can be written as 1/4-1- 

with K > 0. The Fourier series expansion for u{z) is then of the form 

u{z) = 2/^/^ ^ p{n)Kii^{2Tr\n\y)e{nx) 
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with a Bessel function of imaginary argument. We may suppose that our 
cusp forms are eigenfunctions of the Hecke operators T(n) for all positive inte- 
gers n and that u{x + iy) is even or odd as a function of x. Thus T{n)u = t{n)u 
for certain real numbers t{n), which are again called Hecke eigenvalues, and 
u(— z) = ±m(z). Comparing Fourier coefficients on both sides, one may ver- 
ify that p{n) — p{l)t{n) and p{—n) — ±p{n) for all n > 1. The Maass (wave) 
forms Uj constitute an orthonormal set of non-holomorphic cusp forms arranged 
so that the corresponding parameters kj determined by the eigenvalues 1/4: + 
lie in an increasing order. We write Pj{n) and tj{n) for the corresponding Fourier 
coefficients and Hecke eigenvalues. Now let 

= |Pi(l)lVcosh(7rKj), 

and let 

oo 
n=l 

stand for the Hecke L-function attached to Uj. As the Ramanujan-Petersson 
conjecture tj{n) ^ n*^ holds in a mean value sense, the series converges abso- 
lutely for Re s > 1. Again, Hj{s) can be continued to an entire function. We 
have the bound 

E o.,H^{^^^KHoi'>K. (10) 

(See Theorem 3.4.) 

Lastly, in this paper, the following notation will be adopted: Vinogradov's 
relation /(z) ^ g{z) is another notation for f{z) = 0{g{z)). We let e stand 
generally for a small positive number, not necessarily the same at each occur- 
rence. 

2.2 The needed lemmas 

In this section we shall gather some auxiliary results which will be used during 
the course of the proofs of the theorems. 

For estimating the Gamma function we shall frequently use the following 
Stirling's formula without further reference: In any fixed strip h < a < c we 
have 

|r(a + it)\ = y2^|tr-i/'e-l*l"/2(i + 0(|i|-i)), 

for \t\ — >■ oo. For a proof, see Olver [30J p. 294. Further we have 

r(x) = V2^x"-i/2e--(i + 0{x-^)), (11) 

for all positive real x — >■ oo, see [2¥J, Eq. (1.4.25). 

The proof of the following classical estimate for the fourth moment of Rie- 
mann's zeta- function on the critical line can be found e.g. from [33) . p. 147: 

4 

dt-^Tlog^T (12) 
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for T > 1. Along with this estimate we have 



dt < log^^ T 



(13) 



for T > 1 by Heath-Brown [S]. Moreover it is known that 

C,{l + it) 

as \t\ > 1. For a proof, see e.g. [33J, p. 132. 

We have an important tool arising from spectral theory: 

Lemma 2.1 (The spectral large sieve). For K>\,\<lS.<K,M> \ and 
any complex numbers Om we have 



E 



arntj{m) 



2<M 



< {KA + M){KMY ^ \a„ 



i<M 



For a proof, see Theorem 1.1 in [18j or Theorem 3.3 in ^29j. The continuous 
analogy is the following: 

Lemma 2.2. For K real, A > 0, A/ > 1 md any complex numbers a™ we have 

2 



K+A 
K 



arncr2irim)m' 



m<M 



dr < (A^ + M)^^' 



m<M 



uniformly in K, as a2ir{f) — X^dl/*^^""- 
Proof can be found in [31 J, Lemma 1.12. 

A crucial role is played by a spectral decomposition of the shifted convolution 
sum over the divisor function: 

Lemma 2.3. Let f be a positive integer and W a smooth function of compact 
support on (0 , oo) . Then 



^ d{l)d{l + f)W j = -2 m{x; f)W{x) dx+^— 



jl/2 roo 



/-'V2„(/) 



|C(l + 2ir)|2 

oo i?(fc) 

;e 

fe=6 j=l 



e(r;W^)dr + /i/2^a,t,(/)i?| Qj 9(^^,5 PF) 
lf'^'f:T.(~^)'^^^'^t,Mf)Hl, (I) E (k~l;w) (14) 
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Here m{x; f) is defined by (1-12) in f27\/ . cr^ir^f) is as in the previous lemma 
and 

e(.;Ty)4i?e{(l + -^)s(..;I4^)} 

+ + (15) 

with 

xW{x)dx. (16) 

Here F{. ,.;.;.) is the hypergeometric function and {r i— >■ — r) stands for an 
expression similar to the preceding one, but with r replaced by ~r. 

For the proof, see [27J, Theorem 3. 

Finally we introduce the basic inequality in the proof of the classical large 
sieve: 

Lemma 2.4 (Sobolev). Let a < u < a+A and let the function f be continuously 
differentiable on this interval. Then 




uniformly. 

For a proof, see Montgomery [26], Lemma 1.1 applied to 

3 Proofs of the results 
3.1 Theorem [TH] 

We start by adding a set of x log(l/C/) real- valued smooth weight functions 
9^ {.t) Z-sum so that their sum produces an approximation of the char- 

acteristic function of the interval [1,2] with an error of size ^ U and their 
supports widen step by step by factors 2 when we move away from the end 
points 1 and 2. Thus if the supports of the first and last weight functions are 
of length U — L~^+^, then the next ones are of length x 2U and so on, and the 
slopes of the weight functions cancel out each other. Hence we let gs{x) stand 
for a compactly supported function on some interval of length x 6 contained in 
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[1, 2]. Moreover, g{x) = 1 on an interval of length x S and ^'■'^•'(a;) S ^ for 
each 1/ > and a; € M. Here < 5 < 1. Now 



5 < E E E + ^)'^(" + ^ + 

fr^Fn^N 1 = 1 



n + I 



"~2 I iTi{x; f)W{x)dx 



+ N 



l+e 



«^'E E |eiK/;i) + e2(n,/;i) + e3K/;i)|' + A^i+^F, 

f^F n^N 

with W{x) = gs (^^^) ■ Here e^(n, /; L) [v = 1, 2, 3) stands for the {v + l)th 
term on the right hand side of (|14l) . 

The contribution of the term e2{n, /; L) turns out to be the most difhcult to 
bound. Let us denote 



^2 



with e(r; W) as in ^ 
integral in 



EE 

f^Fn^N 



1 

ir, — 
' 2 



As the support of W tends to infinity, we have in the 



ir;l + 2ir; 

X 



1 + VTT: 



-l-2ir 



X F 




(17) 



by the quadratic transformation formula for the hypergeometric function from 
Lebedev [24|, Eq. (9.6.12). We write the RHS of dTTj) by the definition of the 
hypergeometric series. Now clearly the fcth term in the series is non-oscillating 
and of order x"^'^. Therefore it suffices to study only the leading term 1, others 
behaving similarly until some sufficient constant, and the contribution of the 
rest being negligibly small. (See also Convention 2 in [2lJ-) 

By integrating repeatedly by parts we notice that Hj can be truncated to be 
<C N^~^^S~^L~^ with negligibly small error. We further decompose the range 
[1,N^+''S-'^L-'^] into dyadic intervals Kj K with 1 < A' < N^+'^S-'^L-^, the 
number of these intervals being <^ N"^ . Next we prepare ourselves to use the 
duality principle (see for example pp. 169-170 in [H]). Let 6/_„ stand for any 
finite sequence of complex numbers such that J2fr^F Sn-Af = ^- Then 



E E^/^"^^(/)^('*^'^) 



-1 
3 
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with 



zL + n 



-l/2-tK, / 1 + ./l + 



/ 



-1— 2ZKn 



zL-\-n 



gs{z)dz. 



Note that Kj replaced by —Kj can be taken care of with complex conjugation. 
By (dni) this is further 



(18) 



We want to next apply the spectral large sieve, so we shall use Sobolev's 
lemma as done in Jutila's paper [14J on p. 454 to get rid of the parameter 
Kj in the integral I{Kj;W). Therefore the range {K,2K] for kj is split up into 
segments of length A in such a way that the integral I{y; W) remains essentially 
stationary as y runs over a segment. That is 



A log 



V 



4/ 



< logiV. 



In this way, the second term in the upper bound in Lemma 12.41 will be compa- 
rable to the first. An appropriate choice would be A = 1. 

Now we divide the K^-sum in (jlSp into subsums of length A, and apply 
Lemma 12.41 to each subsum arriving at the bound 



N'KL^F-'Y] 

1=0 Ji<+i 



E 



K+l<Kj<K+l+l 



E E bf,nW)ny;W) 



E E^/-*.(/)|:^(y;w^) 

f^F n~N ^ 



dy. 



We next apply the spectral large sieve to the subsum over Kj, and finally 
add the results corresponding to all subsums. This leads us to the bound 

N-KL^F-'{K + F)J2 E^/^"^(2^;^) 



+ 



dy. 



Now we use Cauchy's inequality for the integrals I{y;W) and -§-I{y]W), and 



attach a suitable real- valued smooth weight function u{y) to the y-integral. We 
set u to be compactly supported on the interval [K/2, ZK], u{y) = 1 on [K, 2K] 
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and u^'^\y) <Ci/ K~'^ for each v >0 and y € K. By opening the squares we end 
up with a double sum X^th^at J2n2~N- repeated partial integration over the 
y-integral we now see that we may truncate ni — n2 ^ N^~^^ . Hence we 
finally have an upper bound 

f^Fn~N In-n'I^JVi+^if-i 

and we conclude with the desired upper bound. 

The contribution of ei(n, /; L) to S can be estimated by exactly the same 
argumentation as above. We need the estimation 




following readily of and (|13p by the Cauchy's inequality. Hence 

Si < N^+^LF. 

We could use also a more straightforward argumentation, but these steps make 
it easier to follow the proof of Theorem ll.4l 

Using ([7]) , ([8]) , ([9]) and pTj) the trivial estimates suffice for the estimation of 
e3{n,f;L). 

3.2 Corollary O 

Using straightforward calculation it is easily seen that the difference of our main 
term and that of Atkinson's is 

g /.(2L+n)// 

/ m(x; f)dx - {M{2L + n, /) - M{L + n, /)) < FLN-^+', 

^ J{L+n)/f 

whence the result follows. 

3.3 Theorem Q 

Analogous to the proof of Theorem ll.il For the size of the Fourier coefficients see 
([8]). Instead of Lemma [2731 we use Lemma 4 in [22], originally due to Motohashi 
[28j. The counterpart for pU)) is Lemma 4 in [T^, and instead of ([T^ we need 
Lemma 3 in |14| . 
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